Photovoltage Dynamics of the Hydroxylated Si(lll) Surface Investigated by Ultrafast 

Electron Diffraction 
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We present a novel method to measure transient photovoltage at nanointerfaces using ultrafast 
electron diffraction. In particular, we report our results on the photoinduced electronic excitations 
and their ensuing relaxations in a hydroxyl-terminated silicon surface, a standard substrate for 
fabricating molecular electronics interfaces. The transient surface voltage is determined by observ- 
ing Coulomb refraction changes induced by the modified space-charge barrier within a selectively 
probed volume by femtosecond electron pulses. The results are in agreement with ultrafast photoe- 
mission studies of surface state charging, suggesting a charge relaxation mechanism closely coupled 
to the carrier dynamics near the surface that can be described by a drift-diffusion model. This 
study demonstrates a newly implemented ultrafast diffraction method for investigating interfacial 
processes, with both charge and structure resolution. 
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The problem of directly converting solar energy into 
electrical energy through photocarrier generation has re- 
cently gained tremendous attention due to new types of 
photovoltaic (PV) cells utilizing nanoparticles and molec- 
ular interfacesiii2i£'^£ An important parameter determin- 
ing the efficiency of the solar cell is the survival rate of 
photogenerated carriers reaching the electrodes, which is 
strongly affected by surface recombination and the trans- 
port characteristics at the contacts. In a Schottky-type 
PV junction, the carrier dynamics are driven by a space- 
charge barrier layer £ A photovoltage is generated as the 
photoinduced electron-hole pairs near the semiconductor 
interface move to reduce the band-bending. More gen- 
erally, diffusion of non-equilibrium carriers charging 
of the surface states^ and photoionizationli are differ- 
ent contributors for emfs at semiconductor junctions, at- 
tributing to the overall magnitude and sign of the gen- 
erated photovoltage. For a molecular interface, the exis- 
tence of discrete surface and/or charge transfer states 
are crucial in directing the carrier separation in new 
generation PV cells, such as dye-sensitized TiC>2 meso- 
porous films 2 - and quantum dot solar cellst^ In these 
cases, the filling of a surface state can strongly affect 
the charge distribution in the space-charge layer, modi- 
fying the transport characteristics, as revealed in photoe- 
mission studies i 10 i 12 ' 13 Surface charging is also essential 
for explaining surface photochemical processes, includ- 
ing catalysis^ and molecular transport^ Here, we use 
a novel diffractive potentiometry approach to determine 
the transient photovoltage based on the ultrafast elec- 
tron diffraction (UED) technique, whose ability to re- 
solve surface structure evolution has been demonstrated 
previouslyj 16 ' 17 

We study the photoinduced surface potential changes 
on a Si(lll) surface terminated with hydrophilic hy- 
droxyl (OH) groups, a prototypical system for fabricat- 
ing molecular electronic devicesi 18 ' 19 The experiment is 
performed in an optical pump-diffraction probe arrange- 
ment, in which the femtosecond laser pulse is used to ini- 
tiate carrier generation near the surface and in turn, the 



formation of a transient surface voltage (TSV), while a 
charge-sensitive electron pulse probes the ensuing carrier 
relaxation dynamics. A Coulomb refraction formalism is 
deduced to treat ultrafast diffraction data to determine 
the TSV. Using a laser fluence (F) of 22 mJ/cm 2 , we 
found a band-flattening TSV of 300 mV generated 30 
ps after the laser excitation. At higher fluences (F > 72 
mJ/cm 2 ), the TSV continued to rise, even beyond the 
bandgap energy. This surprisingly large TSV is at- 
tributed to a modified barrier due to a non-equilibrium 
surface charge migration, induced by photoexcitation. 13 
We observed picosecond charge injection and relaxation 
dynamics with timescales similar to results obtained us- 
ing ultrafast photoemission studies of surface carrier dy- 
namics in a vacuum cleaved Si(lll)-(2 x I) surfaced 
This rapid recovery can be understood based on a drift- 
diffusion model that couples surface carrier dynamics 
with transport in the space-charge layer. This suggests 
that the hydroxylated silicon surface posseses a lower 
concentration of trap states, ideal for mediating molecu- 
lar electronic transport study^ 

The hydroxylated Si(lll) surface was prepared ex situ 
using a wet chemistry method. The as received Si(lll) 
wafer (Silicon Quest Intl., p-type, roughness < 12 /im) 
was treated with a modified RCA procedure^ The na- 
tive oxide layer, which usually contains a high density of 
trap states, was removed and an ultrathin (< 5 nm) ther- 
mal oxidation layer with high electric quality was grown. 
First, the wafer was immersed in H2SO4/H2O2 (7:3) so- 
lution at 90 °C for 10 min to remove surface contami- 
nants. The surface was then etched with a saturated 
NH4F solution 21 to remove the native oxide layer, fol- 
lowed by a bath in NH 4 OH/H 2 02/H 2 (1:1:5) at 80 °C 
for 20 min to further remove inorganic residues. Finally, 
a thin layer of oxide that affords a hydroxylated surface 
was grown by immersing the wafer in HCI/H2O2/H2O 
(1:1:6) at 80 °C for 10 min. The wafer was rinsed with 
deionized water (17.8 Mfi-cm) for 10 min at the end of 
each cycle. This procedure yields a density of w 10 15 
cm~ 2 silanol groups (Si-OH) on the surfaced 



2 




FIG. 1: Electron scattering geometry. The incident beam, Si 
(energy eVo, incidence angle 9i, Ewald sphere origin O) scat- 
ters off the Si(lll). Constructive interference occurs when 
the scattered beam (s D ) crosses an intersection of the Ewald 
sphere with a reciprocal lattice rod. Here, surface steps from 
the Si(lll) cleave, tilt the reciprocal lattice rods by cut , re- 
sulting in interference maxima below the projected Zeroth 
Order Laue Zone. From observable peaks (0,0,24), (0,0,27), 
and (0,0,30), we measured As = 1.985 A, implying interlayer 
spacing dm = 2-k / As = 3.165 A and 8 cut ~ 1.35°. The 
ABC stacking of Si(lll) diamond lattice has lattice constant 
9.406 A, implying an interlayer spacing of dm = 3.135 A, in 
agreement with our measurements. 
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FIG. 2: (a) Ray traces of the electron path with and without 
a surface potential (V s ) present. V s induces a shift <5 in the 
diffraction maxima, (b) The shift in diffraction peak maxima 
as a function of O . Lower order diffraction maxima shift more 
than higher orders in the presence of a surface potential. 



The experiment was conducted in the ultrafast electron 
crystallography setup^ employing the reflection geome- 
try to gain surface sensitivity. As described in Fig. [I] the 
incident probe electron has a kinetic energy Kq — eVb, 
where Vq — 30 kV is the accelerating voltage, correspond- 
ing to an incident reciprocal wavevector Si = 90 A" 1 . At 
an incidence angle, 9i = 6.8 °, the Ewald sphere intersects 
the (0,0) reciprocal lattice rod of the Si(lll) substrate 
near momentum transfer wavevector s w 18 A -1 , pro- 
ducing an in-phase (0,0,27) diffraction maximum on the 
screen. Since the Si(lll) substrate has a slight miscut, 
(&cut ~ 1.35°), the diffraction maxima appear slightly 
below the projected Zeroth Order Laue Zone. Weak 
(0,0,24) and (0,0,30)22 diffraction maxima arising from 
the surface roughness were also observed. These maxima 
are associated with the interference pattern of Si(lll) 
planes, evident from As = 1.985 A -1 between neighbor- 
ing peaks satisfying the reciprocal relationship dm = 
2tt/ As with an accepted dm value of 3.135 A. The elec- 
tron penetration depth in silicon is estimated to be rj 2 
nm based on the Scherrer criterion^ using the width of 
the (0,0,27) diffraction maximum. 

To determine the surface potential, we examine the 
Coulomb refraction introduced by a surface potential for 
electron diffraction, as shown in Fig. HJa). The surface 
potential, V s , alters the kinetic energy of the electrons 
submerged in the crystal, K' = e(Vo+V s ), without chang- 
ing the transverse momentum, thereby causing a refrac- 
tion at the interface. The index of refraction can be de- 
scribed by n e — yj (V s + Vq)/V ~ For surface diffrac- 
tion, this alters the electron incidence angle towards the 



lattice plane, thus introducing a net angular shift (5) of 
the diffraction maximum. We calculate the shift arising 
from a surface potential to be 

V, (9 5 - 9,9 Q + 5 2 /2f - g|(go + 8f 

v Q & + e y ' { ' 

where all symbols are defined in Fig. [Ha). In the opti- 
cal pump-diffraction probe experiment, V s is replaced by 
the transient surface voltage, AV S , to evaluate the pho- 
toinduced refraction shift of diffraction maxima at 9 . 
The dependence of S on 9 Q is shown in Fig. [2fb), calcu- 
lated using Eqn.fl]) with our incidence angle (Qi — 6.8°) 
and AV S = 2.7 volts. This refraction shift has a 'non- 
structural' characteristic, exhibiting a more significant 
deviation at smaller values of 9 Q . This differs fundamen- 
tally from structure-induced shifts, which behave oppo- 
sitely. For small angle diffraction the Bragg law yields a 
monotonic increase in 5 as a function of o for a given 
structural change Ad, described by 5/(9 + 9i) = Ad/d, 
i.e. the higher order diffraction maxima experience pro- 
portionally larger shifts than the lower order ones. This 
clear distinction allows us to separate Coulomb-induced 
effects from structural changes by analyzing multiple 
diffraction maxima. To demonstrate the principle, we 
will follow both (0,0,24) and (0,0,27) diffraction maxima 
to examine the TSV dynamics. 

The characteristic features associated with the 
Coulomb effect were observed experimentally as shown 
in Fig. [3Ja). Using a laser fluence of 72 mJ/cm 2 , the 
time-dependent shifts of diffraction maxima (0,0,24) and 
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(0,0,27) were followed by adjusting the arrival time of 
the electron probe pulse relative to that of the exciting 
laser pulse. The shift of the (0,0,24) diffraction maxi- 
mum in reciprocal space is markedly greater than that 
associated with (0,0,27), pointing to a Coulomb effect. 
We expect a negligible structure- induced shift because of 
the small absorption coefficient in silicon (a -1 w 1/xm) 
for the near infrared excitation (hv =1.55 eV), indicat- 
ing that the thermal energy deposition is spread out 
over a large volume. The maximum temperature rise 
in Si can be estimated from ATj, = C^ l aF(\ — R){hv — 
E g )/hv + C^Arie-hEg, where R = 0.37 is the reflectiv- 
ity, Cl = 2.08 x 10 6 Jm~ 3 K _1 is the lattice heat capac- 
ity, and An e -h is the carrier density drop due to Auger 
recombination.— The laser fluences we applied ranged 
from 22 to 72 mJ/cm 2 , which give rise to only moder- 
ate increases in lattice temperature, about 20 to 170 K, 
respectively, above the initial temperature, To = 300 K. 
Hence the changes induced by lattice heating, As/s, are 
on the order of 10 based on the thermal expansion co- 
efficient for Sii2£ In contrast, the carrier generation is 
more significant, ranging from 5 x 10 20 to 1.4 x 10 21 
cm -3 , which corresponds to an increase of several or- 
ders of magnitude from the intrinsic carrier concentra- 
tion (10 15 - 10 16 cm" 3 ) within the first ps^ It is ex- 
pected that such a large increase in carrier density will 
easily flatten out the initial surface band-bending, which 
is w 0.3 V, 27 leading to a change As/s » 7.8 x 10~ 3 for 
the (0,0,27) maximum. Indeed, a band-flattening TSV of 
0.38 V was determined at the lowest fluence applied (F 
— 22 mJ/cm 2 ). However, TSVs higher than the initial 
band-bending were also observed; 1.74 and 2.97 V for 
fluences of 46 and 72 mJ/cm 2 , respectively. In Fig. [3Jb), 
the TSVs deduced from Eqn.fl]) are consistent for both 
the (0,0,24) and (0,0,27) maxima, for all three fluences, 
demonstrating the robustness of the diffractive TSV mea- 
surment. 

The significant photovoltage induced by the high flu- 
ence femtosecond laser cannot be explained merely by 
a barrier-layer effect, where photoinduced e-h pairs near 
the semiconductor interface separate to screen the surface 
charges responsible for the band-bending. Without mod- 
ifying surface charge population, this theory would pre- 
dict a maximum TSV corresponding to the dark current 
band-bending 28 and the decay time of the TSV would be 
comparable to the characteristic bulk recombination time 
(ks 10 — 100 ns for Si). This is not what is observed here. 
In contrast, we have witnessed a much more rapid pro- 
cess, with the TSV surging within 30 ps following laser 
excitation and a decay on the order of 100 ps for each 
of the three different fluences, without any saturation, 
as shown in Fig. [3Jb). Enhanced surface band-bending 
was also observed by Marsi et al. in studying the transient 
charge distribution at the SiC^/Si interface following UV 
free electron laser excitation using time-resolved pump- 
probe core photoemission spectroscopy^ The cause of 
the enhancement is attributed to electron diffusion into 
the thermally grown Si02 overlayer and subsequent accu- 
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FIG. 3: (Color online) (a) Transient shift in diffraction max- 
ima induced by the TSV. The shift in peak (0,0,24) nearly 
doubles that of (0,0,27), indicating the presence of a TSV. 
(b) The TSV for the laser fluences employed. Insets: semi- 
log TSV, demonstrating the contrasting behavior in the TSV 
rise for the given peaks. 



mulation at its surface. The recombination of the excess 
electrons was found to be slower with thicker oxide layers. 
For 12 A oxide thickness, the surface charge recombina- 
tion was found to be comparable to the typical excess 
carrier recombination time in the Si(lll) space-charge 
layer, of the order of 100 ns. For very thin (a few A) ox- 
ide layers, a very effective 'surface' recombination process 
was observed within their time resolution of « 200 ps. 
These observations are consistent with a picture in which 
the surface recombination rate is determined by the sep- 
aration between the excess surface electrons and excess 
holes in the Si space-charge layer. Indeed, much shorter 
recombination time is found for the decay of the surface 
state population in a vacuum-cleaved Si surfaced This 
speed-up in surface recombination for an ultrathin oxide 
layer (< 1 nm) observed by Marsi et al. is attributed to 
the enhanced overlap of bulk evanescent states of Si with 
the surface states. In that case and in the studies re- 
ported here, the electrons accumulated at the thin oxide 
surface are so close to the Si space-charge layer that their 
role in the recombination process is similar to the normal 
surface state. 
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FIG. 4: (Color online) Transient surface state population 
measured by photoemissioni^ compared with the measured 
TSV (rescaled; F — 72 mj/cm 2 ), demonstrating very similar 
decay rates. 



The decay behavior of the TSV on the hydoxylated 
Si(lll) surface was compared with the surface state dy- 
namics on a vacuum-cleaved Si(lll)-(2 x 1) surface, re- 
vealed by a photoemission study (Fig.[4]),m, In this figure, 
the TSV is rescaled to match with the data extracted 
from Ref. 10. We find a striking similarity between the 
two, suggesting that the surface state population is in 
quasi-equilibrium with the bulk carrier dynamics in the 
space-charge layer, so long as there is strong coupling 
between the surface state and the bulk evanescent state. 
Further evidence of the strong role played by carrier dy- 
namics on the surface recombination is revealed by the 
power-law decay of the TSV, shown in Fig. 03a) , where 
an exponent near -1 (-0.93±0.03) is ascribed to the char- 
acteristic decay of surface charges, with a time constant 
of 100 ps. 

To describe the population decay of the surface state 
driven by a space-charge drift recombination, we intro- 
duce a simple two-slab model, illustrated in Fig. [5Jb)- 
(d). The electrons accumulated at the surface lead to a 
rapid rise in surface potential within the thin oxide layer 
that decays slowly into the Si space-charge region, as de- 
picted in Fig. [5jb) . The non-equilibrium space-charge 
layers, with initial separated charge density <7rj, are mod- 
eled as two separate slabs, separated by a distance l e -h 
in a dielectric medium (e w 13 for Si), shown in Fig.[5]Jc). 
The power-law characteristics are manifested in the sim- 
ple rate equation da(t)/dt = a(t)/r r , with the space- 
charge recombination time r r = l e -h/fJ-E depending on 
the transient field E, which is directly related to the sur- 
face charge a(t). Its solution a{t) = ao/(t/r c + 1) has a 
characteristic timescale r c — l e ~h£/ V^Oi which is the time 
for the TSV to drop by 50% from its initial value, and 
corresponds to the induction period in the log-log plot of 
the solution, as described in Fig. Q5]^d) . For t >> t c , the 
t _1 behavior emerges, similar to what is observed exper- 
imentally [Fig. [5Ja)]. To compare with the experimental 
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FIG. 5: (Color online) (a) Log-log plot of the TSV. (b) Sur- 
face potential (14) at the SiOH/Si(lll) interface, depicted by 
solid gray line, (c) Two-slab model employed (schematic) to 
describe the space-charge region, (d) Carrier drift causes the 
surface charge (cro) to diminish with characteristic time r c , 
after which, ao falls as t , obtained from two-slab model. 



results, we assume a linear decay of the surface poten- 
tial barrier (V s ) over the space-charge layer l e -h, and 
determine the size of the induced space-charge regime 
to be l e -h ~ 400 nm (using carrier mobility, /i=100 
cm 2 V -1 s -1 according to Ref. 25, our measured r c of 
100 ps, and V s of 2 V), which is in agreement with the 
laser penetration depth in Si. 

Other mechanisms that would also affect the TSV de- 
cay include the ambipolar diffusion of the bulk excess 
carriers away from the surface and the decay of the sur- 
face state charges into trap states. Here, we can ex- 
clude ambipolar diffusion as the main cause for TSV de- 
cay for two reasons. First, the characteristic diffusion 
time, Tdiff, is significantly longer than what we observed 
( T diff = le^h/D e -h ~ 700 ps, based on the excess e-h re- 
gion generated by the laser, l e -h ~ 1 (tm, and ambipolar 
diffusivity £> e -/i ~ 15 cm 2 s _1 )S Second, a t -1 / 2 decay 
is expected in the case of one-dimensional diffusion, 29 
which is appropriate for photoexcitation when the laser 
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spot size (600 /im) is large compared to the penetra- 
tion depth (rj 1 /im). The rapid decay of the TSV af- 
ter 100 ps precludes any long-lived trap state from play- 
ing a significant role. Halas and Bokor invoked coupled 
drift-diffusion-Poisson equations for bulk carriers to self- 
consistently treat the surface state dynamics.^ Such a 
calculation is able to describe the general behavior of 
the surface state depopulation and TSV decay (to be re- 
ported elsewhere^). However, the essential power-law 
behavior can be elucidated with a simple space-charge 
recombination model, as described above. 

While the simple Coulomb refraction formalism de- 
scribed in Eqn.((T]) seems to be adequate for deducing 
consistent TSVs from diffraction maxima, especially for 
the recombination dynamics, some departure between 
the TSVs deduced for (0,0,24) and (0,0,27) is evident at 
short times, as depicted in the insets of Fig. [3jb). This 
discrepancy is believed to be caused by a non-equilibrium 
barrier layer, which deviates significantly from a homo- 
geneous potential model assumed in Eqn.fl]) within the 
electron probed volume. More rigorous treatment of the 
Coulomb refraction formalism with a non-homogeneous 



potential is needed to fully account for the short-time 
dynamics. 

In conclusion, we have demonstrated a new method 
of directly measuring transient surface potential us- 
ing ultrafast electron diffraction. In comparison with 
other commonly used surface imaging techniques, such 
as Kelvin Force and Scanning Tunneling Microscopy, we 
are able to examine surface potential changes from the 
femtosecond to nanosecond timescale. This technique 
is complementary to ultrafast photoemission methods 
in which the energetics of the photoemitted electrons, 
rather than the potential, are directly measured. Since 
our employed electron probe energy is reasonably high, 
it allows surface potential characterization at high laser 
excitation fluences, which are usually inaccessible by the 
low energy photoemission techniques due to the space- 
charge effect.— We have revealed a regime of strong pho- 
toinduccd surface voltage rise, which can serve as an im- 
pulsive current source for investigating time-dependent 
transport at nanointerfaces. 
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